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CNJ ABSTRACT 

o 

^""^ We observed the Sh 2-233IR (S233IR) region with better sensitivity in near-infrared 

• • than previous studies for this region. By applying statistical subtraction of the back- 

.^H ground stars, we identified member sources and derived the age and mass of three 

X distinguishable sub-groups in this region: Sh 2-233IR NE, Sh 2-233IR SW, and the 

C^ "distributed stars" over the whole cloud. Star formation may be occurring sequentially 

with a relatively small age difference (~ 0.2 — 0.3 Myrs) between subclusters. We found 
that the slopes for initial mass function (F ~ —0.5) of two subclusters are flatter than 
that of Salpeter, which suggests that more massive stars were preferentially formed in 
those clusters compared to other Galactic star- forming regions. These subclusters may 
not result from the overall collapse of the whole cloud, but have formed by triggering 
before the previous star formation activities disturbed the natal molecular cloud. Addi- 
tionally, high star formation efficiency (^40%) of the subclusters may also suggest that 
stars form very efficiently in the center of NE. 



Subject headings: infrared: stars - stars: formation - stars: pre-main sequence - stars: 
individual(Sh 2-233IR) - ISM: jets and outflows - ISM: cloud 



1. Introduction 



Within 2 kpc, 70 — 90% of stars have been found to form in the clusters (Lada &: Lada 2003) 



These newly formed stars, however, are embedded deeply in dense molecular cloud and there exists 
various difficulties in identifying sources and deriving physical parameters. Several statistical tools, 
such as the K-band luminosity function (KLF), the color-color or color- magnitude diagrams, etc., 
have been successfully used to constrain the characteristics of deeply embedded stellar clusters 



(e.g., Lada & Lada 2003). Development of new sensitive equipments, such as, large mosaic IR 
array, enables us to observe embedded stellar cluster in greater detail than previous effort. The 
Sh 2-233IR region is a well studied star-forming region in our Galaxy (e.g., Porras et al.||2006 ) and 
it provides an ideal laboratory for understanding properties of embedded clusters. Therefore, we 
revisit this region with higher sensitivity near-infrared data taken toward a wider region to better 
constrain the properties of the embedded stellar cluster than previous studies. 

Toward the direction of the Galactic Anticenter, Sh 2-233IR (hereafter S233IR), as a part of the 



Sh 2-235 GMC complex (Reipurth Sz Yan 2008), is located at a distance of about 1.8 kpc (Porras 



et al.| 2000), in association with four extended HII regions, Sh 2-231, 232, 233, and 235 (Heyer 
et al.||l996t . Its position coincides with an IRAS source, IRAS 05358-^3543 (02000 = 05''39™10^ 
(^2000 = +35°45'19"). S233IR is classified as a massive star formation region ( Sridharan et al.|2002 ). 



showing CO outflows (Snell et al. 1990) and various maser emissions associated with this region 



(Henning et al. 1992 Tofani et al. 1995 Beuther et al. 2002b; Menten 1991 Minier et al. 2005). 



The K' band image of this region shows many bright stellar sources and also extended nebulous 



features associated with dust emission (Hodapp 1994). The two embedded young clusters, Sh 2- 



233IR SW (hereafter SW, located in the south-west direction from the center) and Sh 2-233IR NE 
(hereafter NE, in the north-east direction), are notable in this region, with remarkable H2 bow 
shocks associated with the NE cluster (Porras et al. 2000). In addition, numerous studies have 



been reported, especially for the NE cluster, including polarimetric observations (Jiang et al. 2001 



Yao et al.|[2000l ), molecular outflows ( [Beuther et "ar]|2007[ [Mao fc Zeng|r2004t [Beuther et al.|[2002a 



Cesaroni et al. 1999 Larionov et al. 1999), and mid-infrared sources (Longmore et al.]|2006) 



In this paper, we revisit this S233IR region with wider field of view, higher resolution, and better 
sensitivity data in near-infrared and radio wavelengths. Our goal is to understand star formation 
history in terms of age, star formation efficiency and initial mass function. We summarize our 
observations at radio and near-infrared in § [2| and the observed results in § [3} Discussion is given 
in § |4j and the summary in § [5} This is the first paper in the series of our work for the S233IR 
and associated region. Here we focus on properties of the embedded stellar clusters in the S233IR 



region. The H2 shock features and cloud kinematics have been studied by Ginsburg et al. (2009), 



and the extended nebulous X-band features, young stellar object and star formation in a larger 



area will be discussed in the second paper (Yan et al. 2010) 



2. Observations and Data Analysis 
2.1. CFHT observations 



Near infrared data were obtained using the Wide-field Infrared Camera (WIRCam, Pascal et al. 
2004) equipped on Canada- France-Hawaii Telescope^ (CFHT) on November 18, 19 and December 
20, 2005 for the H2 and Kg bands, and February 4, 2007 for the J, H, and K continuum bands. Its 
field of view is 20' x 20' with a pixel scale of 0.3" per pixel, and the seeing was 0.5" — 0.7" during the 
observations. Exposure times on each frame were 6 and 30 (58 for H2) seconds for the wide band 
and the narrow band filters, respectively. Total integration times for each filter were 420 seconds 
for J and H, 235 seconds for Kg-, 1755 seconds for H2, and 2400 seconds for Br7 and K continuum. 
The CFHT WIRCam standard pipeline (the IDL Interpretor of the WIRCam Images, aka i'iwi) 
was used for basic image processing, such as, the standard bias subtraction, flat-fielding and bad 
pixel masking. Sky background of each image was subtracted using the comoving-averaged frame 
taken during the observations. Astrometry and photometry corrections were performed against the 
2MASS catalog in order to correct the image distortion and flux level. And then the images were 
analyzed using IDL, IDL Astronomy Library and Phot Vis, a GUI interpretation of IDL DAOPHOT 
package. The completeness and limiting magnitude were estimated by adding artiflcial stars in the 
observed images. The limiting magnitudes were estimated to be 21.0, 20.0 and 19.0 for the J, 
H, and Kg bands, respectively, at a 90% completeness level. The saturation level of each fllter is 
44000 ADUs, corresponding to 13.4 (J), 13.5 {H) and 12.8 [Kg) magnitude, and the photometry 
for saturated stars were adopted from the 2MASS catalog directly. The error of photometry at 
faint end is ~ 0.6 magnitude. Within the observed region, we identifled 567, 715 and 714 stars in 
the J, H, and Kg bands, respectively. 



2.2. ARO observations 

The CO hue ( J = 3 - 2, 345.7959899 GHz) observations were carried out with the SMT0 10 
m telescope at Mount Graham, Arizona, using the dual-channel SIS 0.8 mm receivers operated in 



^The Canada-France-Hawaii Telescope (CFHT) is operated by the National Research Council of Canada, the 
Institut National des Sciences de I'Univers of the Centre National de la Recherche Scientifique of France, and the 
University of Hawaii. 

■^The Submillimeter Telescope is operated by the Arizona Radio Observatory (ARO), Steward Observatory, Uni- 
versity of Arizona, with partial support from the Research Corporation. 



single-sideband dual polarization mode with image sideband rejection of at least 16 dB on April 6, 
8, 10, 25 and 28, 2008. The spectrometers used were a 2048 channel acousto-optical spectrometer 
(AOS) with a spectral resolution of 500 kHz per channel and 1024 channel Forbes Filterbanks 
(FFBs) with a spectral resolution of 1 MHz per channel. All spectrometers were used simultane- 
ously. On-the-fly (OTF) mapping, centered at a = 05^39™10'' 6 = +35°45'19", was used to map 
a 20' X 20' area. The total number of scanning rows was 120 with 10" spacing and a scanning 
rate of 10"/sec. The calibration was done using the absolute position switching mode after fin- 
ishing every individual row. System temperatures varied from 500 to 1500 K during observations 
depending on the weather conditions. The observed data were calibrated to the antenna tempera- 
ture scale, Ta, which was corrected for atmospheric attenuation, and then converted to main-beam 
temperature by T^b = zr^, where the main-beam efficiency, 77^6, is ~0.65 (see the ARO web-site. 



http://aro.as.arizona.edu/, for the details). 



3. Results 

3.1. The S233IR molecular cloud 

The CO J = 3 — 2 emission shows that the molecular cloud associated with S233IR is located 
on the ridge connecting two other star forming regions in the northwest and southeast directions 
about 3.5 pc away from the S233IR cloud at a distance of 1.8 kpc (Fig.[l|a)). The northwest region 
coincides with the HH region Sh 2-233 (IRAS 05351+3549, |Casoh et"ar]|1986[ ), which is excited by 
a B1.5II star (Hunter & Massey 1990). High-velocity CO J = 1 — emission was found toward 



this region (Jiang et al. 2000| >. Towards the southeast of S233IR, there is another star forming 



core called G173.58+2.45 (IRAS 05361+3539, Shepherd k Churchwell 1996), where a smaU YSO 



cluster (Varricatt et al. 2005) and a molecular outflow exists in the east-west direction (Shepherd 



& Watson 2002). 



The main goal of the CO J = 3 — 2 observations was to determine the boundary of the molecu- 
lar cloud embedding the S233IR cluster. Fig. ^a) shows that the dense region of the S233IR cloud 
is centered on the density peak with a roundish shape. The two extended Ks emission features 
near the center of the CO J = 3 — 2 emission correspond to the NE and SW clusters. At the 20 K 
km s~^ level, there is a common envelope for three star formation regions. The average size of the 
30 K km s~^ contour level of the CO ( J = 3 — 2) emission is about 1.5 pc and the half-power width 
(HPW) is ~0.8 pc. We set the boundary of the S233IR cloud with this 30 K km s"^ (r ~2.5') 
contour level and checked the member sources of the embedded cluster within this boundary. 



The column density of molecular hydrogen can be calculated from the velocity integrated 
main-beam temperature (Tmb = -^^)- 

Vmb 

H2 87rz.3 

jnh2 



1 



i Om^E^/kTe: 
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where A^i is Einstein coefficient, Q(T) is the partition function, ^ is the abundance ratio of 
molecular hydrogen and CO, gu is the statistical weight of the states and J (Tex) is defined as 



J{Te.. 
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Eq. ([T]) can be simplified as: 
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where v is in the unit of GHz, A^i = A, 
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33. 2K. Here we use ^ ~ 10^ as 



CO 



suggested by Lequeux (2005) and T^x is assumed to be 20 K (Ginsburg et al. 2009). The peak 
intensity was found at a = 05'*39"'12'' 5 = +35°45'46" with Ta = 107K, and derived column 
density is Nh2 ~ 1.1 x 10^^ cm~^. The total mass of the cloud was estimated to be ~ lOOOMp,. 



IQ. 



We also estimated the column density using the relation derived by Prerking et al. (1982) for 
dust extinction from the average extinction value, Ay ~ 13.0 magnitude (see §3.2), of the S233IR 
cloud, 

Nh2 = Av X 0.94 X lO^^cm"^. (4) 

The derived column density of S233IR is 1.2 x 10^^ cm~^, which gives a total gas mass 1100 Mq, 
within 30 K km s^^ contour (1.5 pc) in diameter. We derived similar values in two methods 
mentioned above, but adopted the result from the dust extinction since the infrared result may 
better correlated with dust emission. 



3.2. Stellar photometry and extinction correction 

S233IR is a complicated region, which contains stellar clusters embedded, extended nebulous 
emissions and shocked outflow gas, as shown in Fig. [lib). Properties of outflow gas in H2 band is 



discussed by Ginsburg et al. (2009) and the extended emissions will be discussed in detail in the 



forthcoming paper. Here we focus on the properties of the embedded stellar sources. The identified 
point sources within the cloud boundary were categorized into 3 groups, the SW cluster, the NE 
cluster and the "distributed sources" associated with the molecular cloud S233IR except of the 
NE (Fig. [l|c)) and SW clusters (Fig. [ijd)). The centers of NE and SW are at a = 05^39""12'' 
S = +35°45'59" and a = 05'^39™09'' 6 = +35°45'10", respectively. The radii of both clusters are 
set to be about 30" based on the morphology and distribution of extended nebulae. 

Fig. ^ shows the luminosity histogram of the identified point sources (solid line) in different 
filter bands (J, H and Ks) with 0.5 magnitude bin. Comparing with previous data by|Porras et al 



(2000), our observations were made with deeper limiting magnitudes (~ 2 magnitude) and better 
image resolution (~ 3 time better), which results in much large number of detected sources. For 
example, a bright source (IR 93) in the NE cluster, which was considered as one point source in 
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Fig. 1. — (a) The CO J = 3 — 2 integrated intensity map in the velocity range of visr = —30 — 
km s~^ {contours), overlaid on the the WIRCam Ks band image of the size 25' x 25'. The outer 
white box is the region observed with the SMT using the on-the-fly mapping mode and the inner 
red box represents the S233IR region. Contour levels of the CO map were shown from 20 to 100 
K km s~^ with spacings of 10 K km s~^. The peak flux density is 177 K km s~^ at the center, 
(b) The color composite image near S233IR made using J (blue), H (green) and Ks (red) images. 
The H2 emission is overlaid with contours on the composite color image. The lowest H2 contour is 
shown as flux level to be 0.27x10-^ mJy. The image center is at a = 05'*39'"10'' 6 = +35°45'19" 
with the size of about 5' by 5'. (c) The color image of NE cluster made using J (blue), H (green) 
and Ks (red) filters, (d) The color image of SW cluster made using J (blue), H (green) and Ks 
(red) filters. 
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Fig. 2. — Histograms of star counts in different filters. The histograms with sohd and dashed 



Unes are our data and data by Porras et al. (2000), respectively. Vertical lines represent the 90% 



completeness limits of this work (dotted line) and of Porras et al. (2000) (dashed line). 



Porras et al. (2000) data, was resolved into several fainter sources in our data. The two histograms 
are almost identical for the sources brighter than ~ 15 magnitude. 

Fig. Is] is the color-color diagram {H — Ks vs. J — H) for the identified sources in the S233IR 
cloud. The stars detected in the SW and NE cluster regions are marked with squares and triangles, 
respectively, and the other sources ("distributed sources") are marked with circles. Filled symbols 
are for the sources within the classical T-Tau stars (CTTSs) range. The fraction of IR excess 
(non-J band detection) stars, which include the stars in the CTTS region or not detected in the 
J band, are 35% (12/34) and 36% (13/26) for the SW and NE clusters, respectively. The CTTS 
and infrared-excess {H — Ks > 1.5) sources are identified in Fig. Hi with red crosses and green 
boxes, respectively. Some very red objects located at outskirts of this molecular seem to be isolated 
star-forming sources. The nature of these sources will be discussed in the second paper. 

We identified about 800 sources in total within the cloud boundary. It was expected that a 
significant number of background sources exist toward the target cloud in the observed bands. It is 
an important but very difficult job to separate the member sources from the background. Without 
proper motion data for the cluster, statistical subtraction is the only way to determine the member 
stars. To subtract the background stars, we chose a region outside the molecular boundary, centered 
at a = 05''38'"40.1'' 6 = +35°42'07.7", as a reference field with no extinction. The reference field 
was chosen in the area well outside the cloud boundary based on radio observations. By assuming 
that there exists the same background stars in the target cloud (the S233IR cloud), we removed 



expected background sources statistically, following the procedure summarized by Jose et al. (2008): 
select a star on the color-magnitude diagram (CMD) of the reference field, eliminate the star 
with same color and magnitude on the CMD of the target field within the observed uncertainties 
{A{H — Kg) < 0.08 and IS.Ks < 0.07). By repeating this procedure for all sources in CMD, we 
subtract background sources in the target field. There were 254 sources cleaned in a 5' radius 
region. This leads to a number of ~ 2 background sources in 30" radius region. In other words, 
92% and 94% of the sources are associated with SW and NE cluster, respective. 

Fig. ^ shows the color-magnitude diagram {Ks vs. H — Ks) for the sources in the S233IR 
cloud, the reference field, and the background cleaned "member sources", respectively. The red 
lines in the panels represent the age-averaged isochrone between 1.0 to 10.0 Myrs from D'Antona 



& Mazzitelli (1994). Although there should exist various uncertainties in this process, we think 
that this method is the most efficient way to determine the member sources at least within the 
observational uncertainties. Additionally, stars on the left side of the isochrone in Fig. [5] were 
removed as foreground stars. The final result, after the background elimination, is shown in the 
right panel of Fig. [5] 

With these "member" sources, we made extinction correction for pre-main sequence stars 
(PMSs) and CTTSs with two different ways, respectively, using two diagrams (Fig. ^ and Fig. ^. 
Extinction corrections for PMSs were made by projecting them back to the mean isochrone along 
the reddening vector on Fig.pl (^H — K g)"^^^'^'"'^'^ = (// — K5')*™'^ + 0.063 x Ay^ by assuming a normal 
reddening law ( Rieke fc Lebofsky|1985 ). The error of dereddened Ks magnitudes is expected < 0.2 



from the mean isochrone we used (Massi et al.||2006). The extinction can be derived by comparing 



observed and intrinsic color of each star. On the other hand, extinctions of the CTTSs (filled 
symbols in the figures) were corrected by projecting them back to the CTTS locus on Fig. pj The 
averaged visual extinction {Ay) were found to be 9.8±5.2 (the SW cluster, 36 stars), 28.9ibl0.4 (the 
NE cluster, 25 stars), and 13.0ibl0.4 ("distributed stars", 124 stars). These extinction are larger 
than the previous values by a factor of about 1.3 — 2 (PCS), which result from deeper observations 
and also from more reliable analysis by using both diagrams (Figs.pl and pi) in deriving extinction 
amounts. 



3.3. KLF modeling 

For decades, the most reliable method in optical astronomy used to determine the age of a 
cluster is using the HR diagram (HRD), comparing the positions of member stars with theoretical 
evolutionary tracks on the HRD. For embedded young cluster, this method is not validated because 
most of the members are not in the main sequence stage and can be observed only in longer 
wavelengths than optical. The K-band luminosity function (KLF) is a simple tool to study the 



properties and estimate the age of an embedded cluster (Lada k, Lada 2003 Yasui et al. 2006) 



Definition of tlie KLF can be expressed by the following equation: 



dN 



dN 



dlogM^ 



druk dlogM^: dmu 



(5) 



where nik is the K-band luminosity and M^, is the stellar mass ( Lada &: Lada|[2C)03 > . The first term 
in the right hand side is the underlying stellar mass function and the second term is the mass- 
luminosity relation (MLR). It is noticed that the KLF of clusters peak at different magnitude, 



depending on the difference between their ages and star formation history (Muench et al. 2000). 



Simple Monte Carlo simulations were carried out to construct the model KLFs. The simulation 
was done in three step. The first step is to assume an initial mass function (IMF). Two IMFs were 



used in our simulation, they are Trapezium IMF (Muench et al. 2000) and the IMF from Miller 



& Scalo (1979, hereafter MS79). Then, we convert the mass function to the luminosity function 



using a mass-to-luminosity relation (MLR) from the isochrones of the PMS models(D'Antona & 



Mazzitelli 1994 1997, 1998: Siess et al. 2000). The stellar luminosities were finally converted to 



the K-band luminosity, mx with bolometric correction (Flower 1996) and stellar intrinsic color 



correction (Bessell & Brett 1988). By repeating this procedure with different age inputs, we fit 



model KLFs to observed KLFs and estimate the ages of clusters. In Fig.|6j we show a result derived 
for embedded clusters in S233IR, tp be discussed in the following section. 



3.4. Age estimates from KLFs and derived IMFs of the embedded clusters 

Ages of the clusters embedded in the S233IR cloud were estimated by comparing the observed 



and model KLFs (Yasui et al. 2006). Using the method mentioned in ^3.3, we constructed model 



KLFs for ages from 0.07 to 50 Myrs. By comparing the model and observed KLFs (Fig. pi), we 
derived the ages of ~ 0.5±0.1, ~ 0.25±0.1, and ~ 1.5±0.3 Myrs for the SW, NE, and the distributed 
stars, respectively. Although the MLR for the age younger than 1 Myr is uncertain between different 
PMS models ( Baraffe et al.|[2002 Yasui et al.||2006p , this method allows us estimating cluster ages 



quantitatively. The age difference of these stellar groups was first noticed by Porras et al. (2000) 



based on the J-band data analysis. It is understood that there are uncertainties in estimating 
cluster ages because we detect 25 and 36 IR sources in NE and SW cluster, respectively. However, 
the model KLFs are not sensitive to the initial mass functions (IMFs) used. Therefore, the current 
results do suggest that relative age differences exist among these clusters roughly with the amounts 
mentioned here. Using the estimated ages of the clusters, the luminosity function was converted to 
the mass with the model isochrones ( [D'Antona k Mazziteilil|1994| |1997| |1998[ [Siess et al.||2000[ ) by 
applying the bolometric correction ( Flower 1 1996 ) to the expected stellar intrinsic colors (Bessell &; 
Brett||1988D. The total masses of >46, >30, and ~110 Mq were derived for the NE, SW, and the 



distributed stars, respectively. 

With the determined age, the IMFs of the embedded clusters can be calculated using the 
mass-luminosity relation and its slope F, which is defined as dN/dlog{m) = M^ (cf. Muench et al. 



10 



2002), to constrain the cluster properties. The size of the IMF bin was set to be 0.5 Mq, which 
is larger than the uncertainty of mass estimation propagated from luminosity error. We derived 

r 0.42 ± 0.11, ~ -0.10 ± 0.08, and 0.46 ± 0.07 for the SW , NE, and distributed stars, 

respectively. Although the derived F values for the NE and SW clusters are uncertain due to small 
numbers (25 source for NE and 36 sources for SW), it still indicates that the F values in this region 
are much flatter than the Salpeter slope. This may suggest that more massive stars have efficiently 
formed in the S233IR region compared to other star forming regions in our Galaxy, such as the 



sources listed by Muench et al. (2002); Figueredo et al. (2002, 2005); Leistra et al. (2005, 2006); 



Jose et al. (2008); Harayama et al. (2008); Pandey et al. (2008), etc., who have derived slopes 



F ~ — 1 in young embedded stellar clusters. 



3.5. Star formation efficiency 

The star formation efficiency (SEE) is defined as the ratio of the total stellar mass to the total 
stellar and gas mass (cf. Wilking & Lada 1983). Typically, the SEE ranges from ~10% to 30% 



( Lada Sz Lada|[2003 see Table 2 and reference there in). SEEs have been found low in young or low 
mass star formation regions (eg. Serpens, Rho Oph and NGC 1333) but increased to about 30% 



in more evolved massive star- forming regions (Lada & Lada 2003). We have determined the total 
gas mass of S233IR, ~ 1100 Mq, and the embedded cluster mass, ~ 180 Mq (§.[3lJ). The SEE is 
~17% for the S233IR region. Our analysis covers the whole S233IR region and the member stars 



were determined more carefully than previous studies for this region (for example, Beuther et al. 
2002a|bt [Kfao k Zeng||2004l ) . 



It should be noticed that the derived SEE may represent the lower limit, since significant 
number of sources have not been identified because of the extended emission features and high 
extinction environments in the central parts of the embedded clusters, this caveat also applies to 
F. The latest star-forming cores are found in the NE cluster. Young massive stars are still deeply 



embedded in the dense molecular cores (Beuther et al. 2007; Leurini et al. 2007) and associated 



with maser emissions and energetic shocked H2 gas (Menten 1991; Porras et al. 2000 Ginsburg et 



al. 2009). Because of the age differences among stellar groups, it is possible that star formation 



has not resulted from the overall collapse of the S233IR molecular could at the same time. If stars 
form continuously, the single SEE value for the whole cloud may not represent a measure of the 
final gas-to-stars efficiency. We therefore investigated SEEs at different radii in the S233IR cloud. 
Since the total gas mass in local areas is difficult to separate from the whole cloud, the SEEs were 
derived by comparing the stellar mass and the total gas mass along the ring, centered at the cloud 
center, of the projected width of about Ar = 0.2' (projected distance of ~0.1 pc), as shown in 
Eig.JT} The center is at a = 05''39'"11.8'* 6 = +35°45'51.1", which corresponds to the peak of the 
integrated CO emission. As expected, the SEE changes with the location in the cloud. High SFE 
values up to > 40% are found in the central region (r < 1'). This radius covers the region of the 
NE and SW clusters and both of them represent later stages of star formation. The high value may 
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suggest that there are local enhancements in massive star formation, possibly triggered by nearby 
star formations within the S233IR cloud. 



4. Discussion 



Based on the radio and extinction map ( Heyer et al.||1996 Reipurth fc YanpOOS ) , the S233IR 
molecular cloud is a part of the much larger Auriga GMC complex, which is associated with several 
large HII regions and active star-forming regions. The large scale (~ a few degrees) morphology of 
this GMC complex looks very clumpy and is composed of several HII regions, such as Sh 2-231, 232, 



233 and 235 (Reipurth & Yan 2008). The CO emission shows a filamentary and arc- like structure 
associated with three HII regions. At the northeast direction, an HII region, Sh 2-231, is ~ 10' 



away from S233IR (Mao Sz Zeng 2004). In addition, there are many newly formed OB stars, such 



as those in Aur OBI, associated with this GMC complex (Reipurth Sz Yan 2008, Figure 1). The 



S233IR clump appears to be the latest of star formation events. 

Three groups of stars were identified in S233IR region, SW, NE and the distributed stars. 
The background contamination in distributed stars were eliminated by a statistical method. The 



age differences among these subgroups were first noticed by Porras et al. (2000) using the J-band 
analysis. Based on better data quality of our work, we derived ages ~0.5 Myrs for SW, ~0.3 
Myrs for NE, and ~1.5 Myrs for the distributed stars. The distributed stars could be the first 
generation of stars formed in this cloud. The SW and NE clusters, with the relative distance 
of about 0.5 pc at a 1.8 kpc distance, are later generations. The age differences between stellar 
groups suggests sequential star formation, which has also been found in other star-forming regions. 



such as NGC 7538 (Ojha et al. 2004), Sh 2-247 (Puga et al. 2008), Sh 2-157 (Chen et al. 2009) 



and NGC 3576(Purcell et al. 2009). It is interesting that the S233IR cloud, a typical star-forming 



cloud in size (r ~1.5 pc) and mass (~ 1100 Mq), contains sub-clusters formed in different epochs 
and star formation still continues without disrupting its natal molecular cloud. Stars appear to 
be forming continuously in the youngest NE cluster. For the distributed stars, the derived IMF 
slope, r ~ —0.5, is flatter than in other star-forming regions in our Galaxy, which indicates that 
more massive stars have been formed preferentially in this cloud than in other typical star-forming 
regions. The F value of the NE cluster, the latest star-forming sub-core in this cloud, is even flatter 
than the distributed stars, suggesting that the probability of massive star formation is higher at 
this region. This result is consistent with the radio observation toward the center of NE cluster 



(Beuther et al. 2002a). Eventually the cluster will disperse the natal molecular cloud and appear 



as visible stars. Once the massive stars are formed, the radiation from massive star will clear the 
molecular could and stop the star formation process in the core region. 

The enhancement of the local SFE is also found in the central dense core as we derived in 



^3.5 although we have to take various uncertainties into consideration, for example, the uncertain 
total H2 density, the cloud boundary, member star determination, undetected sources especially 
near the central part of the sub-clusters showing saturated infrared emission, etc. The local high 
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SFE ( ^40%) suggests that the later generation of stars can be formed very efficiently in dense 
cores in the same molecular cloud. In the outskirt of the cloud, where the first stellar generation 
is located, the SFE is consistent with a typical molecular cloud. The trigged star formation may 
have enhanced the formation of massive stars, which results in a high SFE and flatter IMF slope 
in the central region. It is known that the NE cluster is the latest star forming core in S233IR, 



which is associated with the powerful energetic features, such as shocked H2 gas flows (Ginsburg 



et al. 2009). These energetic activities are thought to be a part of massive star formation process. 
The S233IR region provides an interesting example of active and energetic massive star formation, 
which could be further studied with next generation instruments. 



5. Summary 

We revisited the S233IR (Sh 2-233IR) region with wider field of view and more sensitive near- 
infrared data than previous studies. There exist three distinguishable sub-clusters in this region: 
the NE and SW clusters, and the distributed stars. Using the color-color (H — Ks vs. J — H) and 
the color-magnitude {Ks vs. H — Ks) diagrams, we subtracted background stars statistically and 
identified 25, 33, and 151 stellar sources for the NE cluster, the SW cluster, and the distributed 
stars, respectively. We derived ages ~0.3, ~0.5, and ~1.5 Myrs, and mass 45, 30, and 107 M©, for 
SW, NE and the distributed stars, respectively. The flatter IMF slope (F > —0.5) for subclusters 
suggests that more massive stars were preferentially formed in S233IR. The first stellar generation 
is the distributed stars located in the GMC complex in which S233IR is embedded. Afterward, the 
sequential star formation in this cloud may have been triggered before its natal molecular cloud 
is disturbed by previous star formation events. The averaged SFE of this region is ~17%, and 
>40% at the center of S233IR, the NE and SW clusters. This result suggests that stars form very 
efficiently in local dense cores within the cloud center. 
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Fig. 3. — The H — Kg vs. J — H diagram of the identified point sources in the S233IR region. 
Sources found in the SW, NE cluster regions, and the "distributed sources" are marked with squares, 
diamonds, and circles, respectively. CTTS are shown with filled symbols in the gray-shaded region. 
The blue dashed line is the CTTS locus and the solid blue lines the track of main sequence stars 



(Bessell & Brett 1988 Meyer 1996). Red dashed lines indicated the direction of extinction, and 



the red arrow is for A 
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Fig. 4. — Location of the sources identified as CTTSs (red crosses) and infrared-excess sources 
(stars detected in H and Ks, green boxes). The background color image is the same as Fig. fltb). 
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Fig. 5. — The color-magnitude {Ks vs. H — Ks) diagrams. Source symbols are the same as in 
Fig. [3] Red lines represent the age-averaged isochrone between 1.0 to 10.0 Myrs from D'Antona &: 



Mazzitelh (1994) 
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Fig. 6. — Histograms of K-band luminosity function in the SW, NE, and the distributed stars as 
indicated in the panels. The solid and dash-dotted lines represents the model KLFs based on the 



IMF of Miller & Scab (1979) and Muench et al. (2002), respectively. 
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Fig. 7. — The radial dependence of star formation efficiency in S233IR region. The SFEs were 
derived by comparing the stellar mass to the total gas mass along the ring, centered at a = 
05''39"'11.8* 6 = +35°45'51.1", with annuh of width Ar = 0.2'. 



